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Abstract. This paper examines line profile evolution due to the linear expansion of circumstellar 
material obsverved during a microlensing event. This work extends our previous papers on emission 
line profile evolution from radial and azimuthal flow during point mass lens events and fold caustic 
crossings. Both "flavours" of microlensing were shown to provide effective diagnostics of bulk motion 
in circumstellar envelopes. In this work a different genre of flow is studied, namely linear homologous 
expansion, for both point mass lenses and fold caustic crossings. Linear expansion is of particular 
relevance to the effects of microlensing on supernovae at cosmological distances. We derive line profiles 
and equivalent widths for the illustrative cases of pure resonance and pure recombination lines, modelled 
under the Sobolev approximation. The efficacy of microlensing as a diagnostic probe of the stellar 
environs is demonstrated and discussed. 
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1. Introduction 

Recently gravitational microlensing has been shown 
to be a powerful tool for probing not only the nature 
and distribution of dark matter but also the astro- 
physical properties of the source being lensed (see, 
for example, Gould 2001 for a comprehensive review). 
While the vast majority of microlensing events ob- 
served to date can be adequately approximated as 
point sources lensed by a point mass lens, there is a 
small but significant fraction of events in which the 
source must be modelled as extended, which has im- 
portant consequences for the evolution of the broad- 
band lightcurve and spectrum during the microlens- 
ing event and produces observational signatures that 
are now readily detectable - chiefly due to the intro- 
duction of microlensing 'alert' networks. This devel- 
opment has allowed intensive photometric and spec- 
troscopic monitoring of many microlensing events, 
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spearheaded by the PLANET collaboration (Albrow 
et al. 1998). The alert strategy is especially well- 
suited to fold caustic crossings - in e.g. close binary 
lensing events - since observation of the first caus- 
tic crossing permits prediction of the second caustic 
crossing, allowing intensive follow-up observations to 
be scheduled. Moreover, every fold caustic crossing 
must be treated as an extended source event, making 
them ideal probes of the astrophysics of the source. 

In the past few years 'alert response' broad-band pho- 
tometric observations of extended source microlens- 
ing events have shown clear detections of stellar limb 
darkening on the photospheres of a number of stars, 
including: K giants (Albrow et al. 1999; Albrow et 
al. 2000; Fields et al 2003), a metal-poor A dwarf 
in the Small Magellanic Cloud (Afonso et al. 2000), 
a G/K subgiant (Albrow et al. 2001a) and an F8- 
G2 solar-type star (Abe et al. 2003). These observa- 
tions have been used to estimate limb darkening pa- 
rameters and to constrain stellar atmosphere models. 
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comparing for example the LTE ATLAS models of 
Kurucz (1994) with the PHOENIX 'Next Geireration' 
models of Hauschildt et al. (1999). 
In addition Albrow et al. (2001b) used the VLT to 
carry out observations of Ha equivalent width varia- 
tion in the atmosphere a K3 bulge giant during a fold 
caustic crossing event. These authors showed that 
highly resolved observations of the microlensed spec- 
trum could be a powerful discriminant between differ- 
ent stellar atmosphere models, and indeed could also 
be a useful method for breaking the near-degeneracy 
between the parameters of different lens models. This 
work essentially vindicated the theoretical modelling 
of Heyrovsky, Sasselov and Loeb (2000), which pre- 
dicted that with 8m-class telescopes it should be pos- 
sible to use highly time- and spectrally-resolved ob- 
servations of microlensing events to carry out 'tomog- 
raphy' of stellar atmospheres. 

A notable feature of the theoretical literature on 
extended source microlensing events, however, has 
been the comparative neglect of extended circum- 
stellar envelopes. Although the model developed in 
Coleman l)1998|l and Simmons et al. (2002) includes 
such a scattering envelope, and Coleman et al. (1997) 
also considered the case of a non-spherical envelope 
(e.g., a Be disk), these computations were carried out 
only for the broad-band photometric and polarimet- 
ric response to a microlensing event. 
In two preceding papers, Ignace & Hendry (1999) and 
Bryce, Ignace and Hendry (2003) calculated the mi- 
crolensing signatures of emission line profiles from 
expanding and rotating spherical shells, lensed by a 
point mass lens and a fold caustic respectively. Their 
results were for highly simplistic velocity fields - con- 
stant expansion or constant rotation - in order to iso- 
late the effects of the microlensing and demonstrate 
clearly the diagnostic potential for deriving velocity 
information about the shells. The authors showed 
that, while in the absence of lensing the integrated 
line profiles for expanding and rotating shells yielded 
identical 'flat top' profiles in each case, microlensing 
clearly breaks the degeneracy between these cases. 
In this third paper, we choose to focus on the par- 
ticular case of circumstellar media in homologous ex- 
pansion, with radial velocity v(r) oc r. This case is 
interesting for two reasons. Firstly, this type of veloc- 
ity law applies to the early phases of nova and super- 
nova (SN) explosions. Since the latter are observed to 
cosmological distances, it becomes increasingly likely 
that such events will be lensed as the light from the 
SN makes its way from a distant galaxy to the Earth 
(c.f. Dalai et al. 2003). Thus a consideration of ho- 
mologous expansion is relevant if only to investigate 
definite signatures of microlensing from emission pro- 



file shapes observed in SN events. Indeed Bagherpour 
et al. (2004) consider the microlensed lightcurves of 
Type la SN, and Bagherpour et al. (2005) extend 
their analysis to consider the impact of microlensing 
on P Cygni profiles in the spectra of Type la SN. 
Secondly, in this paper we employ standard Sobolev 
theory to model the emission line profiles; the Sobolev 
approximation is valid for the highly supersonic 
fiows commonly found in SNe and stellar winds. 
Furthermore the case of homologous expansion leads 
to a significant simplification of the expressions and 
a consequent substantial reduction in computational 
efi'ort. Thus it is worthwhile studying this case to bet- 
ter appreciate the characteristics of how microlensing 
will affect emission lines formed in winds with more 
general velocity laws. 

The structure of this paper is as follows. In Sectional 
we describe the basic features of line formation the- 
ory, and the Sobolev approximation which we employ. 
In Section 3 we derive, and briefly discuss, unlensed 
proflles for the case of pure resonance and pure re- 
combination lines, in a stellar wind undergoing ho- 
mogolous expansion. In Section 4 we present illus- 
trative results, calculating the time evolution of line 
proflle and equivalent width for pure resonance and 
pure recombination lines lensed by a point mass lens 
and a fold caustic. Finally in Section 5 we discuss our 
results and their future applicability in a number of 
specific astrophysical contexts. 

2. Line formation theory 

We first need to derive the emission profile for the 
unlensed case, expressed in terms of intensity as a 
function of frequency and position within the enve- 
lope of circumstellar material. For an envelope in bulk 
motion such that the flow speed greatly exceeds the 
thermal broadening, the locus of points contributing 
to the emission at any particular frequency in the 
line profile is confined to an "isovelocity zone" (c.f., 
Mihalas 1978). These zones are determined by the 
Doppler shift formula, namely 

= ^^0 (l - ^) , (1) 

where the observer's coordinates are {x, y, z) with the 
line-of-sight along z, is the Doppler shifted fre- 
quency, and Vz = — v(r) • z is the projection of the 
fiow velocity onto the line-of-sight (as indicated in 
Figures n] and 121) . Note that both Cartesian (a;,y, z) 
and spherical (r, iJ, kp) coordinates can be defined for 
the star. Employing the Sobolev theory for line pro- 
file calculation in moving media reduces the radia- 
tion transfer in a moving envelope to a calculation 
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Fig. 1. A schematic showing the contribution from dif- 
ferent regions to the line profile shape, for an observer 
on the right hand side. Much of the line emission is pro- 
duced around the circular tube indicated in projection by 
the dashed lines. On the far side (left hand side) of the 
tube the emission will be occulted by the star. On the 
near side, the circumstellar material will attenuate the 
continuum emission from the central photosphere. 



where 1} is the angle between the radial direction and 
the line of sight. 

Thus isovelocity zones for constant Vz are seen to be 
plane surfaces. In fact, given the spherical symmetry 
of our model, the isovelocity zones are circular plane 
surfaces of radius 



Pn 



(4) 



as depicted in Figure |3 

Consider a projected element of the stellar wind at 
fixed impact parameter p, and angle a measured 
clockwise from x - i.e. x — pcosa and y — psina. 
The emergent intensity Iu{p)^ due to emission from 
the wind is given by 



K,y{r) p{r) Su e dz' . 



(5) 




where is the opacity, p the density. Si, the source 
function, and t^, the optical depth at radial distance 
r, satisfying = + z'^. 

The Sobolev approximation is to evaluate the optical 
depth only at the point where the Doppler shifted 
frequency of a gas parcel exactly equals the observed 
frequency. This delta function response in frequency 
simplifies considerably the expression for the optical 
depth, namely 



Fig. 2. A schematic showing the isovelocity zones (sur- 
faces of constant v^) for our velocity model. These are 
circular plane surfaces which are oriented normal to the 
observer's line-of-sight. We assume that the photosphere 
is of radius 7?, and the envelope is of maximum radial 
extent rmax. 

in which one need consider only distinct isovelocity 
zones. Equation ^ will therefore prove crucial for 
relating the variable profile shape to the kinematics 
of the envelope. 

For the assumed linear expansion, and for an enve- 
lope of maximum radial extent rmax, the velocity v 
at radial distance r is given by 



(2) 



where f is a unit vector in the radial direction and 
I'max is the wind speed at r = rmax- The velocity 
shift as seen by the observer, situated as in Figure ^ 
is then 



-Wmax(?'/?'max) COS = -Wmax (^^/^'max) , (3) 



n^{r) p[r) 5{vz - Vz')dz' 

K.^{r) p{r) Aq 
\dvz/dz\z 



(6) 



The denominator on the right hand side is the line- 
of-sight velocity gradient. For a spherical flow, this 
gradient is given by the expression 



dvz ^dv ^ V 

dz dr r 



(7) 



For the case of homologous expansion, with v{r) = 
■^max (''/'''max), the linc-of-sight velocity gradient sim- 
plifies to dvz/dz — Wmax/''max, which is a con- 
stant. Therefore, although the Sobolev optical depth 
will generally depend on angle d through the factor 
dvz/dz^ homologous expansion is the single exception 
for which the optical depth depends on radius only. 
The optical depth is now 



K^{r) p{r) Ao 



(8) 



^ Note that the spherical symmetry of our model means 
that the emergent intensity is independent of a. 
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Knowing the optical depth as a function of r, we can 
then obtain an expression for the intensity, [p) , us- 
ing equations (4) and (6) and the equation 



P 



In Sobolev theory the intensity then reduces to the 
form 



(10) 



Allowing for line emission from resonance line scat- 
tering and coUisional de-excitation, and defining e as 
the ratio of the coUisional de-excitation rate to that 
of spontaneous decay, the source function can be de- 
rived to be 



(11) 



The two parameters Pc and P are respectively the 
penetration and escape probabilities, which are given 

by 



/3c 



and 



1 

47r 



(12) 



(13) 



where fi* is the solid angle subtended by the star at 
radius r. 

The optical depth is generally anisotropic, but as al- 
ready noted for homologous expansion the optical 
depth is isotropic in this special case. As a result, 
the penetration and escape probabilities are easily 
calculable, and in particular Pc = W(r) P, where the 
dilution factor 



Wir) = 0.5 (1 - y/l-R'^/r'^) 



(14) 



So for linear expansion, we have the well-known 
result that the Sobolev optical depth is a function of 
radius only. Since the source function also depends 
only on radius (or equivalently on impact parameter, 
p, at fixed Vz), so too does the intensity. To compute 
the total emergent flux at a given Doppler shift in 
the line profile due to emission in the circumstellar 
envelope, we must integrate over emergent intensity 
beams as given by 



1 



ly {p) pdpda, 



(15) 



where D is the distance from the Earth. Substituting 
in for the intensity, the integration for the flux be- 
comes 



(9) F,{vz) 



27r 
Z32 



■,(r) (l-e-^-M) pdp. (16) 



where we have, of course, performed the (trivial) inte- 
gral over the angle a, and the limits Pmin and Pmax are 
discussed below. For homologous expansion we know 
that the Sobolev surfaces are disks oriented trans- 
verse to the line-of-sight (i.e., with z — constant). 
From equation ^ this means that pdp ~ rdr, and 
so the above flux integral could be equivalently re- 
formulated as 



2tt 



St, 



(r) (l e-^-('')) rdr. (17) 



We now include the contribution of continuum ra- 
diation from a pseudo-photosphere at radius R. It 
is instructive to consider separately redshifted and 
blueshifted wavelengths, which correspond to > 
and v{z) < respectively. 

For v{z) > (the redshifted side) the total flux 
Ftot{vz) at frequency v and line-of-sight velocity Vz 
may be written as 



FtotiVz) ^ Fi{Vz) + F2iVz) ^ 



where 
Fi{vz) 
and 



27r 



^phot 



pdp, 



(18) 



(19) 



s. 



,(r) (l-e-^-W) pdp. (20) 



Here /phot is the continuum intensity of the pseudo- 
photosphere. If /phot is a constant then the integrals 
involving /phot become very straightforward. Fi rep- 
resents the flux coming directly from the pseudo- 
photosphere while F2 accounts for photons scattered 
by the surrounding envelope which emerge along the 
line of sight. Note that the lower limit of the F2 in- 
tegral is R since the region for p < /? is occulted by 
the pseudo-photosphere of the star. The upper limit. 



satisfies the relation 



Pn 



(21) 



For v{z) < (the blueshifted side) the total flux con- 
sists of three terms, i.e. 



Ftotivz) = Fsivz) + Fiivz) + F^ivz) , 



(22) 
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where 



3.1. Pure resonance lines 



2n 



Pn, 



Fsiv^) = — / /phot pdp, 



(23) 







and 



/pi- 



Here is defined as 



Plim = 



Vi?2 _ ^2 for < z < i? 







for z> R 



(25) 



(26) 



Note tliat if z > i? then the F3 integral is identi- 
cally zero; in this case the flux contribution from 
the pseudo-photosphere is represented fully by the 
term, which also takes account of attenuation by 
intervening wind material along the line-of-sight. 
Finally we must stress that, in accordance with the 
standard notation adopted in the microlensing liter- 
ature, all distance scales are in fact taken as angular 
distances that arc normalised to the angular Einstein 
radius of the lens. This implies that the "fluxes" 
defined in the above equations have rather unusual 
units, however, the results of our line profile calcu- 
lations in the following sections will be displayed as 
ratios normalised by the continuum, so that the (non- 
standard) units of flux will cancel. 

3. Unlensed profiles from homologous 
expansion 

Before including the effects of gravitational lensing 
we first consider the unlensed profiles derived by ap- 
plying the model of the previous section to the il- 
lustrative cases of a pure resonance line and a pure 
recombination line. We have computed line profiles 
for a range of optical depths, assuming a spherically 
symmetric density distribution for the ejecta given 

by 



p{x) = po x~ 



(27) 



where x = r/R. The factor accounts for mass 
continuity within the shell as material expands. The 
scale parameter po denotes the density at the ra- 
dius of the pseudo-photosphere and specifies the total 
amount of gas ejected in the explosion, via the equa- 
tion 



ejecta 



Aur^ pir) dr . 



(28) 



For the case of a pure resonance line the parameter 
e is set to zero in eq. and the source function 

becomes 



= W{r) 



(29) 



The opacity for resonance line scattering does not 
depend explicitly on density, but an implicit depen- 
dence can arise through, for example, the ionization 
fraction of whatever atomic species is being consid- 
ered. For simplicity, we suppress any radial depen- 
dence of the opacity in our model, and parametrize 
the Sobolev optical depth as 



= To (p/po) = To a; 



-3 



(30) 



An especially convenient parameter used to charac- 
terize different line calculations is the line integrated 
optical depth defined as follows 



T = 



Tydx = iro 



1 



1 



(31) 



Thus, for fixed Xmax, it is trivial to compute the value 
of To - and hence the model profile - corresponding 
to a given value of T. 

FigureElshows some examples of model profiles com- 
puted for T = 0.3,1.0,3.0 and 10.0, with stronger 
lines corresponding to larger values of T. In all cases 
a fixed value of Xmax = 5 was adopted, so that from 
equation T = 0.48ro. Note that our chosen value 
of Xmax is responsible for the change of slope of each 
line profile at Vz = — 0.2winax- This is due to the oc- 
cultation of wind material which lies directly behind 
the pseudo-photosphere of the star (at x — 0.2xinax 
in this case). 

Increasing the value of Xmax would result in qualita- 
tively similar line profiles but with a narrower emis- 
sion peak and absorption trough. We can understand 
this behaviour as follows. From equations lO and © 
we can write 



Pi 2 2 

^ "^max 



(32) 



where w = Wz/fmax- From equation (|30|l it then fol- 
lows that for w ~ the optical depth - and hence the 
line flux - is insensitive to x^ax (The range of inte- 
gration in equation Hlt)|) increases as Xmax increases, 
but this has little effect on the line flux because the 
optical depth is proportional to x~^). 
For w ~ ±1, on the other hand, at fixed p the optical 
depth decreases sharply as Xmax increases - resulting 
in reduced line emission and absorption and thus a 
narrowing of these features in the line profile. 
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3.2. Pure recombination lines 

We have also computed model profiles for the case 
of a pure recombination line. As such, scattering is 
ignored, and the source function is assumed to arise 
from an LTE process, so that 



(33) 



for temperature, T. For simplicity, the envelope shall 
be taken as isothermal. The opacity for a recombina- 
tion line is K oc p, and so the scaling of the optical 
depth will for this case be 



To [p/pof 



(34) 



As a p^ process, the optical depth will go as - a 
strong function of radius. 

It is again useful to introduce the integrated optical 
depth parameter given by 



T = 



To [p/pof dx , 



(35) 



so that again, for fixed Xmax, it is trivial to compute 
the value of tq corresponding to a line of given inte- 
grated optical depth. 

Figure 01 shows profile shapes again computed for 
s^max = 5 and for the same range of values of T 
as in Figure |3| From equation (|35|l it follows that 
T w 0.2ro. Note that in contrast to Figure 01 the 
recombination lines do not exhibit any blueshifted 
absorption. That is because we assume that the con- 
tinuum emission is given by a blackbody with the 
same temperature as the envelope. Consequently for 
sightlines that intercept the continuum forming sur- 
face, one has = B^er'^ + By{l — e^'^) = B^. Hence 
the combined emission from line plus continuum will 
always equal or exceed the continuum level when in- 
tegrated over all sightlines. 

As was the case for pure resonance lines, increasing 
s^max would producc narrower profiles but would have 
negligible effect on the value of the peak emission at 
w = Q for fixed T. 

4. Lensed profiles from homologous expansion 

We now consider the effect of microlensing on emis- 
sion line profiles for a stellar wind in homologous 
expansion. Our expression in equation (|15|l for the 
emergent flux is now replaced by an integral over 
emergent intensity beams weighted by the (dimen- 
sionless) lensing magnification A{d), where d is the 
projected separation between the lens and an area 
element of the wind in the source plane, i.e. 



1 

Z)2 



. T = 0.3, 1, 3, 10 









V /v 



Fig. 3. Model line profiles from pure resonance line scat- 
tering, calculated using the method described in the text, 
for four values of the integrated optical depth, T. The 
stronger lines correspond to larger values of T. 




Fig. 4. Model line profiles from pure recombination line 
scattering, calculated using the method described in the 
text, for four values of the integrated optical depth, T. 
The stronger lines correspond to larger values of T. 



Iu{p)A{d)pdpda . (36) 



Note that, while we continue to assume that I^, de- 
pends only on impact parameter, the lensing magnifi- 
cation is a function of both p and a, so that we must 
perform a double integration to compute the emer- 
gent fiux at each radial velocity. Moreover, since the 
magnification pattern of the lens evolves with time, 
as the relative projected positions of the lens and 
wind change, the computed line profile will also be 
time-dependent. 

The form of the magnification function is determined 
by the nature of the lens. In this paper we consider 
two illustrative cases: those of a point mass lens and 
a fold caustic. 
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4.1. Point mass leasing event 



In the typical microlensing situation of a point mass 
lens the (time-dependent) magnification factor takes 
the well-known form (see e.g. Pacyhski 1986) 



A{u) 



u^/u^ -I- 4 



(37) 



where u is the impact parameter: the projected sep- 
aration between the lens and source element, given 

by 



u{t) = 




(38) 



where t^ is the time of closest approach between lens 
and source element, which occurs at w = uq, the min- 
imum impact parameter, and tE is the characteristic 
timescale of the event. Both u and Uq are expressed 
in units of the angular Einstein radius of the lens, 
which is defined as 



4GM Ds - Dl 



DsDi 



(39) 



where M, Ds, Dl are the lens mass, source distance 
and lens distance respectively. The timescale tE may 
be written as 



ti 



Mrol 



(40) 



where /iroi is the relative proper motion of the lens. 
In the limit of large separations the magnification is, 
of course, equal to unity while for small u the mag- 
nification goes approximately as \/u. Thus if a point 
lens transits an extended source, while formally the 
magnification factor is infinite for the source element 
directly behind the lens, when one integrates over the 
finite area of the source an finite value for the mag- 
nification is obtained. 

Figure El illustrates schematically the geometry of a 
point lens event. The outer edge of the circumstel- 
lar envelope (of radius r^ax) is shown in projection, 
together with the trajectory of the lens. When the 
lens is at position i, an area element of the wind at 
position S (defined by polar coordinates p and a - 
see eq. 36) lies a projected distance u from L. 

4.1.1. Pure resonance lines 

We computed continuum lightcurves, line equivalent 
widths and line profiles for the case of a pure reso- 
nance line lensed by a point mass lens. Illustrative 
results are shown in Figure El for a transit with mini- 
mum impact parameter equal to zero - i.e. assuming 




Fig. 5. A schematic view of the source plane for a point 
lens event, showing in projection (i.e. with the observer 
in the z direction) the outer edge of the circumstellar 
envelope (of radius rmax) and the trajectory of the point 
lens, L. An element of the source, at S, with position 
defined by polar coordinates p and a, is shown; it lies at 
a projected separation u from the lens. 



the lens crosses exactly over the centre of the cir- 
cumstellar envelope^ - and adopting a photospheric 
radius R = 0.03 This value is representative of 
source radii determined for observed microlensing 
events exhibiting finite source effects. Indeed Macho 
event M95-30, which was the first well-documented 
example of a finite source event, was found to have a 
somewhat larger radius of i? w 0.075 9e (Alcock et al. 
1997). In the results shown in Figure Owe truncated 
the radius of the circumstellar envelope at Xmax = 5 
and the line profiles were taken to have integrated 
optical depth T = 3. 

The upper left hand panel of Figure El shows the 
continuum lightcurve of the lensing event - i.e. the 
change in apparent magnitude, computed from the 
ratio of lensed to unlensed continuum flux, as a func- 
tion of time. The time axis is measured in units of t^; 
- the lensing timescale introduced in equation (40). 
As expected, the maximum magnitude change occurs 
at minimum impact parameter, when the lens tran- 
sits the centre of the photosphere, and the continuum 
flux is boosted by a factor of more than 60. Only the 
central portion of the light curve is shown; although 
lensing still produces a continuum flux magnification 
of ~ 5 when {t — to)/tE = ±0.2, the effect of lensing 
on the line profiles is essentially negligible for larger 
impact parameters. 



In view of the symmetry of our model the orientation 
of the lens trajectory during transit is irrelevant 
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Fig. 6. The lightcurve, equivalent widths and line profiles produced by a point mass lens transiting a star with 
photospheric radius R = 0.03 9e surrounded by a circumstellar envelope exhibiting pure resonance line scattering. 
The upper left panel shows the change in apparent magnitude of the lensed continuum flux, compared with the 
unlensed flux, as a function of time. Below it the time evolutioin of the normalised equivalent width of the resonance 
line is shown. The four line profiles in the right hand panel correspond to the epochs represented by dotted vertical 
lines in the right hand panel. Thus, for example, profile 'b' corresponds to the epoch of maximum magnification 
(when the lens is at its minimum impact parameter) and is also when the largest change in equivalent width occurs. 
The solid lines in the right hand panel show the lensed line profiles while the dashed lines show the unlensed profiles, 
and are included for comparison. The offset between the profiles is provided for clarity. All profiles were calculated 
for an integrated optical depth of T = 3.0 and assuming a^max ~ 5. 



The lower left hand panel of Figure El shows the time 
evolution in the normalised equivalent width of the 
resonance line. The equivalent width was calculated 
according to the formula 

EW=/^\-l, (41) 

where the angled brackets denote the mean, or ex- 
pectation, value obtained by integrating over Vz- 



The vertical dashed lines labelled a, b, c and d in 
the upper and lower left hand panels correspond to 
the four epochs in the right hand panel at which the 
lensed (solid curve) and unlensed (dashed curve) line 
profiles are compared. For clarity a constant vertical 
offset is introduced between each of the four epochs; 
note that this does not indicate any change in mag- 
nification. 

We can see from the lower left panel that the un- 
lensed value (i.e. for large impact parameter) of the 
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Fig. 7. The lightcurve, equivalent widths and line profiles produced by a point mass lens transiting a star with 
photospheric radius K = 0.03 Qe surrounded by a circumstellar envelope exhibiting pure recombination line scattering. 
The three panels present results in the same manner as in Figure|S] All profiles were again calculated for an integrated 
optical depth of T = 3.0 and assuming a;niax = 5. 



the normalised equivalent width is slightly negative 
because of the asymmetry in the unlcnsed line pro- 
file: the area of the blueshifted absorption trough is 
slightly larger than that of the redshifted emission 
peak. As the event proceeds and the lens crosses the 
outer edge of the envelope, at {t — t[))/tE = —0.15, 
the normalised equivalent width at first shows a clear 
increase, as the lens differentially magnifies a por- 
tion of the isovelocity zones in emission. This stage 
of the lensing event corresponds to profile (a) in the 
right hand panel. As the projected lens position ap- 
proaches the photosphere, however, the absorption 
component of the line profile quickly begins to dom- 
inate. This is because the lens is differentially mag- 
nifying the attenuated continuum emission coming 



from the photosphere. The equivalent width reaches 
its lowest value when the lens has zero impact param- 
eter - i.e. when it lies directly in front of the centre 
of the photosphere. The line profile at this epoch is 
shown in the right hand panel as profile (b). The 
lensing event thereafter is a mirror image of what is 
observed on the incoming trajectory: the equivalent 
width rises sharply as the lens exits the photosphere, 
reaches a peak value close to the edge of the circum- 
stellar envelope before tending to its unlensed value 
at large impact parameter. The two further line pro- 
files at epochs (c) and (d) illustrate this behaviour. 
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4.1.2. Pure recombination lines 

We also computed continuum lightcurves, equivalent 
widths and line profiles for the case of a pure recom- 
bination line lensed by a point mass lens. Illustrative 
results are shown in Figure [7| for exactly the same 
wind and lens parameters as in Figure El 

The upper left panel of FigureQis, of course, identical 
to that of FigureEl since it shows the magnification of 
the continuum flux. The time evolution of the equiva- 
lent width also shows qualitatively similar behaviour 
to that of the resonance line case, although the mag- 
nitude of the maximum change in equivalent width 
is approximately five times larger; this is essentially 
due to the stronger dependence of the optical depth 
on radius for a recombination line. 

Again we see that the equivalent width first increases 
as the lens transits the circumstellar envelope, due 
to the differential magnification of the emission flux 
from the region around the star. When the lens tran- 
sits the photosphere, on the other hand, the differen- 
tial magnification of the attenuated continuum flux 
from the photosphere again causes a sharp decrease in 
equivalent width, which reaches its minimum value at 
zero impact parameter. The symmetry of the isove- 
locity zones and the lens trajectory then results in 
symmetric behaviour during the second half of the 
event. 



4.2. Fold caustic I en sing event 

For a binary lensing event the magnification pattern 
in the source plane is, in general, a rather compli- 
cated function of the lens masses and separations. As 
we remarked in Section 4.1, for a single point lens the 
magnification A{u) goes as 1/m as u ^ 0, so that the 
magnification is formally infinite at only one point 
- coincident with the position of the point lens it- 
self, for which u = 0. For a binary lens, on the other 
hand, the magnification is formally infinite along one 
or more closed curves, or caustic structures, in the 
plane of the sky. Any source crossing into or out of a 
caustic will experience a high degree of magnification, 
and - no matter how small - must be modelled as an 
extended source; i.e. the total magnification must be 
computed by dividing the source into (formally in- 
finitesimal) area elements, calculating the magnifica- 
tion experienced by each area element, and summing 
(or, formally, integrating) the results. Moreover, as a 
source (or an area element thereof) enters a caustic 
structure, two extra images of it are produced by the 
binary lens. These extra images cause a sudden and 
sharp increase in the total magnification. 



If the angular size of the source is very small com- 
pared with that of the closed caustic curve, then gen- 
erally we may approximate by a straight line the por- 
tion of the caustic curve in the vicinity of the source. 
This situation is referred to as the fold caustic ap- 
proximation, and the magnification function in this 
case can be adequately described as (Schneider et al. 
1992) 



A{d) = 



Ao + ^ (d<0), 



An 



{d > 0). 



(42) 



Here Aq is the total magnification of the 3 images 
which form when the source is outside the caustic 
structure; to a very good approximation this term 
is constant during the caustic crossing. The distance 
d is the projected source-caustic separation (i.e. the 
perpendicular distance between the source and the 
line representing the fold caustic) in units of the an- 
gular Einstein radius of the binary lens (which is 
given by eg. 1391 as for a point lens, but with the mass 
M equal to the combined mass of both binary com- 
ponents). The constant bo depends on the parameters 
of the lens system, but is of order unity for typical 
caustics, and in our calculations for simplicity we set 
it exactly equal to unity. 

Note that the sign of d is important here: the excess 
magnification from the extra two images occurs only 
when the source (or an element thereof) lies inside 
the caustic structure. In equation l|42|) . therefore, the 
magnification is defined according to the sign con- 
vention that d increases from left to right, and the 
caustic interior lies to the left of the fold caustic (i.e. 
for d < 0). This situation is illustrated schematically 
in Figure|Sl which shows in projection the outer edge 
of the circumstellar envelope (of radius rmax) and a 
small portion of a closed caustic structure, approxi- 
mated by a straight line. An element of the source, at 
S, with position defined by polar coordinates p and 
a, is also shown; it lies at a projected perpendicular 
distance d from the fold caustic. Note that the fold 
caustic will not, in general, be aligned with the di- 
rection of motion of the source (assumed in Fig. |Hlto 
be along the a;-axis) so that the caustic may sweep 
across the source at an oblique angle. This oblique 
case is shown in Figure |H1 in the examples presented 
in the next two sections, however, we consider only 
the case where the fold caustic is perpendicular to 
the direction of motion of the source. 
Note further that the fold caustic approximation will 
generally be good, since the size of the caustic struc- 
ture will typically be much larger than the size of the 
source and hence curvature of the caustic will not be 
important. The approximation will break down if the 
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Fig. 8. A schematic view of the source plane for a fold 
caustic crossing event, showing in projection (i.e. with the 
observer in the z direction) the outer edge of the circum- 
stellar envelope (of radius rmax). The bold line represents 
a small portion of a closed caustic structure, approxi- 
mated by a straight line. An element of the source, at 
S, with position defined by polar coordinates p and a, 
is shown; it lies at a projected perpendicular distance d 
from the fold caustic. Note that the fold caustic will not, 
in general, be aligned with the direction of motion of the 
source (assumed here to be along the a;-axis) so that the 
caustic may sweep across the source at an oblique angle. 



caustic crossing occurs close to a 'cusp' in the caus- 
tic structure, but we do not consider that case here. 
We will extend our treatment to more general caustic 
models (and to other expansion laws) in future work. 

4.2.1. Pure resonance lines 

We next computed continuum lightcurves, line equiv- 
alent widths and line profiles for the case of a pure 
resonance line lensed by a fold caustic as the source 
exits the caustic structure. Illustrative results are 
shown in Figure|^for a star with pseudo-photospheric 
radius R — 0.003 and assuming Xmax = 5. The 
model line profile shown has T = 3. 
The upper left hand panel of Figure El shows the 
change in apparent magnitude, computed from the 
ratio of lensed to unlensed continuum flux, as a func- 
tion of time (again measured in units of ig). Note 
that we consider here only the effect of the excess 
magniflcation due to the extra images produced by 
the caustic structure; in other words we are assuming 
that the term Aq in equation (|42(l is a constant during 
the caustic crossing, and can therefore be subtracted 



from the total magnification, leaving only the differ- 
ential effect of the two extra images. Hence the mag- 
nitude change due to the extra images alone drops to 
zero after the photosphere has fully exited the fold 
caustic. 

The lower left hand panel of Figure IHl shows the time 
evolution of the line equivalent width, calculated as 
before. We see that the equivalent width first rises 
gently as the fold caustic begins to cross the circum- 
stellar envelope and differentially magnifies the line 
emission from the wind relative to the continuum. 
The increase in equivalent width is at first modest, 
however, because the absorption components of the 
line profile also receive some differential magnifica- 
tion at this stage. Indeed once the fold caustic be- 
gins to cross the star, the dominant effect of the 
differential magnification is to enhance the attenu- 
ated continuum flux from the pseudo-photosphere 
- thus resulting in a sharp dip in the equivalent 
width which reaches its minimum value at epoch (b) . 
Thereafter the equivalent width begins to rise dra- 
matically. This is because an increasing fraction of 
the pseudo-photosphere now lies outside the caustic 
structure, so that the relative effect of the lens on the 
circumstellar material which still lies inside is signif- 
icantly enhanced. This effect reaches its peak shortly 
after epoch (c), by which time the photosphere lies 
entirely outside the caustic but a portion of the cir- 
cumstellar envelope still lies inside, and thus contin- 
ues to be differentially magnified due to the existence 
of the two extra images. As the fold caustic then 
sweeps across the rest of the envelope, the equiva- 
lent width drops sharply again, although at epoch (d) 
the equivalent width is still somewhat larger than at 
epoch (c) - as is evident in the profiles shown in the 
right hand panel. 

4.2.2. Pure recombination lines 

Finally we computed continuum lightcurves, line 
equivalent widths and line profiles for the case of a 
pure recombination line lensed by a fold caustic as 
the source exits the caustic structure. Illustrative re- 
sults are shown in Figure^]for eaxctly the same lens 
and source parameters as in Figure |51 
The upper left panel of FigurelTUlis again identical to 
that of FigureEl From the lower left panel we see that 
the time evolution of the line equivalent width shows 
qualitatively similar behaviour to that of Figure IHl 
although again the magnitude of the change in equiv- 
alent width is significantly larger for a recombination 
line than for a resonance line ~ as we also found for 
a point lens event. The stronger signature of lensing 
is also clearly seen in the right hand panel of Figure 
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Fig. 9. The lightcurve, equivalent widths and line profiles produced by a fold caustic crossing of a star with photo- 
spheric radius R = 0.003 0b surrounded by a circumstellar envelope exhibiting pure resonance line scattering. The 
three panels present results in the same manner as in Figures El and All profiles were again calculated for an 
integrated optical depth of T = 3.0 and assuming a^max ~ 5. 



IIUI where the lensed line profiles at epochs (a) to (d) 
show significantly greater enhancement of the emis- 
sion near hne centre than the enhanced emission at 
the corresponding epochs of Figure |51 or in Figure 
for a recombination hne lensed by a point mass lens. 

Thus we can see collectively from these illustrative 
figures that the time evolution of the equivalent width 
as the lensing event proceeds can be an effective di- 
agnostic of the radial extent of the circumsteUar en- 
velope. The equivalent width variations will also be 
sensitive to the integrated optical depth of the line; 
as we saw in Figures 01 and 01 increasing T resulted 
in a stronger line. Therefore, by studying the time 
evolution of equivalent width for a number of differ- 
ent lines, microlensing can in principle be a powerful 



probe of the run of optical depth as a fimction of 
radius in the circumstellar envelope. 

5. Discussion 

fn this paper we have modelled the effects of mi- 
crolensing for variations of line profiles formed in cir- 
cumstellar envelopes. Although the microlensing is 
quantitatively accurate, the underlying wind model 
includes a number of simplifying assumptions: the 
flow is assumed to be spherically symmetric; reso- 
nance and recombination lines are both treated, but 
without detailed consideration of ionization or tem- 
perature gradients in the flow; the wind velocity law 
is treated as a linear function of radius; the envelope 
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Fig. 10. Same as for FigureEl but now shown for a pure recombination line. 



is truncated at some maximum radius; Sobolev the- 
ory is employed for the line formation. The intention 
of our modeling has been to assess relative effects 
from the line variations for an illustrative model, one 
that has generic application to stellar winds, but that 
would nonetheless require modification when applied 
to any particular case. 

At this juncture it is useful to comment on the di- 
verse kinds of winds observed across the Hcrtzsprung- 
Russell Diagram (hereafter, "HRD"). Referring to 
the introductory stellar winds book by Lamers & 
Cassinelli (1999), most stellar winds fall into two pri- 
mary classes, these being radiatively driven or pres- 
sure driven. Of the former, line-driven winds are 
the primary consideration for hot OB stars (e.g., 
see review by Kudritzki & Puis 2000) as basically 
described by the CAK theory of Castor, Abbott, 
& Klein (1975). Line-driving is also relevant for 



the central stars of planetary nebulae that show 
Wolf-Rayet-like spectra (Tylenda, Acker, & Stenholm 

1993) . The Main Sequence hot stars exhibit absorp- 
tion line spectra from their photospheres; however, 
strong and weak P Cygni lines are abundant at UV 
wavelengths (Howarth & Prinja 1989; Snow et al. 

1994) . In the evolved hot stars, recombination lines 
appear at IR wavelengths (e.g., Lenorzer et al. 2002). 
Also, a substantial fraction (around a quarter, al- 
though the value is uncertain) of B stars are classi- 
fied as emission-line or "Be" stars (Porter & Rivinius 
2003). These have circumstellar disks that produce 
copious recombination lines in the IR, but more in- 
terestingly show Ha in emission, which would be 
relevant for ground-based follow-up monitoring pro- 
grams. Unfortunately, hot stars are relatively rare 
and so are extremely unlikely to be background 
sources for microlensing events. 
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Working toward cooler temperatures in the HRD, 
the winds from Main Sequence stars A-M are not 
well-studied, in part because their mass loss is so 
much weaker than early- type stars. The yellow hy- 
pergiants have more substantial winds, but these 
stars are extremely rare (de Jager 1998). A sub- 
class of the A stars are strongly magnetic - the Ap 
and Am stars. These represent ^ 10% of the A- 
star class (MacGregor 2005), and are stars with sur- 
face magnetic fields in the multi-kilo-Gauss range, 
often modelled as oblique magnetic rotators in which 
the magnetic field is predominantly dipolar (Stibbs 
1950). The wind diagnostics, particularly the por- 
tion trapped by the global dipolar magnetic field 
loops, are likely to be found in the X-ray band, 
as the trapped wind streams from opposite hemi- 
spheres are forced by the magnetic field into super- 
sonic collisions to generate shock- heated gas (e.g.. 
Babel & Montmerle 1997; Townsend & Owocki 2005). 
Again, these types of stars are not often expected in 
finite-source microlensing events. The solar-type star 
observed in the extended source event reported in 
Abe et al. (2003) presumably has a solar-like stellar 
wind, with a feeble mass-loss rate of order 10~^^ Mq 
yr~^. Consequently, wind spectral features accessi- 
ble to ground-based follow-up would be essentially 
non-existent for such a source; instead, coronal diag- 
nostics would be found in the UV and X-ray bands. 
Although one might expect to see effects in Ca ii H 
and K that are sensitive to chromospheric emission, 
our model line profiles are not appropriate for the 
quasi-static chromosphere. 

The most likely class of stars to be found as sources 
in microlensing transit events, either by a single or 
a binary lens, are the red giant stars, a consequence 
of their being reasonably common and fairly large in 
size. Red giant winds have modest mass-loss rates of 
order 10~^° or 1O~^M0 yr~^; however, even these 
estimates are not certain because the origin of their 
winds is not well-understood. The absence of X-ray 
emission for K and M giants (Linsky & Haisch 1979) 
indicates that their winds are not Parker-like, al- 
though they do show chromospheric features. The 
strong Mg II h and k lines and Fe II lines can be 
used to probe cool star winds; however, these lines 
are found in the UV which is not amenable to ground- 
based follow-up. It may be that a combination of 
molecular opacity and some dust formation drives the 
winds of red giants (Jorgensen & Johnson 1992), and 
perhaps certain molecular lines could be used as wind 
diagnostics. 

Red supergiants (RSG) and asymptotic giant branch 
(AGB) stars have more interesting wind features that 
could be studied in microlensing transits. These are 



stars with high mass-loss rates of order 1O~^M0 yr~^ . 
They are radiatively driven, but unlike the early-type 
stars that are line driven, these stars are continuum 
driven because of substantial dust formation in their 
dense and cool winds. Simmons et al. (2002) have 
explored the possibility that scattering polarization 
could be used to probe their winds, but of course 
this is a continuous opacity. Similar to the red giants, 
UV lines can be used to probe the wind flow. Some of 
RSG and AGB stars also show maser emission lines at 
radio frequencies (Elitzur 1992; Habing 1996; Lewis 
1998). These are thought to form primarily in shells. 
Consequently, the model lines presented in this paper 
are less likely to be relevant, although the consid- 
erations of Ignace & Hendry (1999) might be more 
appropriate, if suitably modified for the emissivity 
function of the maser. Once again, stars of the RSG 
and AGB classes are fairly rare. 
There are two other classes of sources for which our 
line profile models are more relevant. One is the ac- 
tive galactic nuclei (AGN) that produce strong wind 
flows (Arav et al. 1995; Proga, Stone, & Kallman 

2000) . The (rest frame) UV spectra of quasars show 
some quite strong P Cygni profiles, such as in C 
IV, suggesting wind speeds in excess of 10,000 km 
s^^. The flows are likely line-driven similar to early- 
type stars. Already, some authors have modelled 
the characteristic emission line profile effects that 
would be observed for different source models during 
microlensing events (Popovic, Mediavilla, & Muiioz 
2001; Abajas et al. 2002), ranging from Keplerian 
disks to relativistic disks to simple shells and jets. 
The microlensing under consideration was by a single 
deflector; however, our spherically symmetric wind 
models are not directly applicable to these cases that 
mostly involve non-spherical geometries. Some have 
claimed the detection of microlensing effects relevant 
to probing the AGN accretion disk (e.g., Chae et al. 

2001) , although none have reported effects pertinent 
to a wind component. We will not discuss this class 
further, and have mentioned it only because the gen- 
eral consideration of emission profile variations dur- 
ing microlensing events involving wind flows has been 
considered in the context of AGN. 

The second class of objects arc supernovac (and the 
likely related gamma-ray bursters). The explosions 
of stars spew gaseous ejecta into space. There is a 
phase in which the ejecta follows homologous expan- 
sion, not because the flow is accelerating, but because 
the ejecta is moving outward with a distribution of 
speeds. Still, the flow can be modeled as a wind, 
and the spectra at different phases shows numerous 
wind-features, ranging from strong P Cygni lines, to 
recombination lines, to forbidden lines (Filippenko 
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1997). Because supernovae (SNe) are so intrinsically 
luminous at early times, they can be seen to large 
distances, and so the probability of lensing increases 
owing to the considerable cosmological path lengths 
involved. 

Recently Dalai et al. (2003) have discussed the im- 
pact of weak lensing due to large scale structure 
on the use of SNc as standard candles. Other au- 
thors have considered microlensing as a tool for 
studying the cjccta of SNe and even gamma-ray 
bursts (GRBs). For example, Schneider & Wagoner 
(1987) considered polarimetric variations of pseudo- 
photosphorcs produced in SNe during a microlens- 
ing event. Baghcrpour, Branch, & Kantowski (2005) 
modelled the impact of microlensing on P Cygni line 
profiles in the spectra of Type la SNe. Gaudi, Granot, 
& Locb (2001) have invoked microlensing to explain a 
brightening event in the lightcurve of GRB-000301C. 
However, thus far a more comprehensice study that 
systematically explores the effects of microlensing on 
emission line features formed in SN and GRB ejecta 
remains to be made; our illustrative results presented 
here are a further step in that direction. 
With respect to SNe, there are several points of cau- 
tion to be made when applying the results of this pa- 
per. First, our models assume a steady-state source. 
Of course, the luminosities, ejecta optical depth, and 
photospheric radii evolve with time in real SNe. Our 
results are applicable to SNe if the lensing event is 
relatively fast compared to photometric and struc- 
tural changes of the ejecta. Under what conditions 
will this be true? Suppose wc take the duration of 
the lensing event to be one week. If the fastest ejecta 
shell moves at 10^ km s~^, then the extent of the 
ejecta shell would be at least r,„ax ~ 40 AU at early 
phases of the explosion in which homologous expan- 
sion is applicable. Moreover, the evolution must be 
such that the lens moves across the ejecta shell (in 
projection) before the dimensions of the photosphere 
change significantly. This means the proper motion 
of the lens should greatly exceed that of the photo- 
spheric expansion. Given the 10** km s~^ ejecta speed 
as a maximum, and assuming a transverse velocity 
of 10^ km s"'^ for the lensing mass, we require that 
£>SN ^ -Dl- Consequently, if the lens were located 
at kpc distances (e.g., in the halo, the bulge, or the 
Large Magellanic Cloud), then the SN should be lo- 
cated at Mpc distances (e.g. at the distance of the 
Virgo cluster). 

As a consistency check, we compare the angular size 
of the SN ejecta 6sn to 9e, as given by 

' max 

^/{Ri./D^) (1-£»l/£'sn)' 



If we assume the lensing object to be typical of 
A4ACH0S of about O.SAf©, then the Schwarzschild 
radius of the lens i?L will be of order 1 km. If we 
further assume rmax ~ 40 AU, and Z^sn ^ D\,, then 

OsN ^ 100 kpc ^^^^ 

where Dl = 10 kpc was assumed. Thus, SNe at 
distances of several to dozens of Mpc would have 
6sn/0e ~ 0.01 — 0.1 at early phases in the SN when 
homologous expansion remains valid. These happen 
to be values typical of our model simulations. 
A further assumption of our models has been the 
common "core-halo" treatment, in which we can ap- 
proximate the photosphere as a "hard" spherical 
boundary, and the wind features are all formed exte- 
rior to this boundary. Recent radiative transfer sim- 
ulation of Type II SNe by Dessart & Hillier (2005) 
show this not to be true. Instead, the radiative trans- 
fer in SNe ejecta is more akin to that of Wolf-Rayet 
stars, for which the continuum forms in the wind flow. 
The purpose of our models has not been to repre- 
sent rigorously any particular type of stellar wind or 
class of stars, but to illustrate the generic effects that 
could be observed with intensive follow-up programs 
for source transit microlensing events. In so doing, 
models for two common types of lines have been con- 
sidered - resonance P Cygni lines and recombina- 
tion emission lines, for both point lensing and caus- 
tic crossing events from binary lensing. Whether the 
source is a hot star or a cool star, whether the lines 
are in the X-ray band or the radio band, the general 
properties of the line variations in terms of lead/lag 
times relative to the photosphere, of equivalent width 
variations, and changes in line profile shape are qual- 
itatively all to be expected. We imposed the homolo- 
gous expansion as a representative flow velocity law. 
Observationally, the velocity law is something that 
would preferably be determined from the line vari- 
ations as the microlensing event evolved. With the 
photospheric crossing time determined from the pho- 
tometric variations, the variations of the line equiv- 
alent width and profile changes with time could be 
converted to radius, calibrated by the photospheric 
radius, to deduce flow velocity and opacity variations. 
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